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Abstract 

We calculate the form factors V, A\, A2 and Aq appearing in the D — * <j> transition 
in light cone QCD sum rule method. We compare our results on these form factors 
with the current experimental results and existing theoretical calculations. 
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1 Introduction 



Semileptonic decays of mesons containing charm and beauty quarks constitute a very im- 
portant class of decays for studying strong and weak interactions. These decay modes of 
heavy flavored mesons are much more clear samples compared to that of the hadronic decay 
modes, since leptons do not participate in the strong interaction. 

Therefore, the study of these decays is one efficient way for determining the elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, as well as for understanding the origin 
of CP violation which is related to the structure of the CKM matrix in Standard Model 
(SM). 

An accurate determination of CKM matrix elements, obviously, depends crucially on 
the possibility of controlling the effects of the strong interactions. For exclusive decays, 
where initial and final states of hadrons are known, the main job is to calculate various 
transition form factors, which involve all the long distance QCD dynamics. So, some non- 
perturbative approach for estimating the long distance effects is needed. Several methods 
have been used to treat these effects, such as quark model, QCD sum rules, lattice theory, 
chiral perturbation theory, etc. Among these approaches, QCD sum rules occupies a special 
place, since it is based on the very first principles of QCD. 

The method of QCD sum rules [1] has been successfully applied to wide variety of 
problems of hadron physics (see [2, 3] and references therein). In this method, physical 
observables of hadrons are related to QCD vacuum via a few condensates. The semileptonic 
decay D — > K°ev e was firstly studied in QCD sum rules with 3-point correlation function 
in [4]. This method, then, successfully extended to study other semileptonic decay decays 
of D and B mesons, i.e., D + — > K°e + u e , D + — > K°*e + u e [5], D — > 7rez/ e , D — ► pev e [6], 
B -> D(D*)£u e [7] and B -> n£u e [8]. 

However, this method inherits some problems, the main one being that some of the form 
factors have nasty behavior in the heavy quark limit itlq — > oo. In order to overcome the 
problems of the traditional QCD sum rules, an alternative method, namely light cone QCD 
sum rules (LCQSR) was developed in [9] and is regarded as an efficient tool in studying 
exclusive processes which involve the emission of a light particle. 

The LCQSR is based on the operator product product expansion (OPE) near the light 
cone x 2 ~ 0, which is an expansion over the twist of the operators, rather than the dimen- 
sions as in the traditional QCD sum rules. All non-perturbative dynamics is parametrized 
by the so-called light cone wave functions, instead of the vacuum condensates in the tradi- 
tional sum rules, which represents the matrix elements of the nonlocal operators between 
the vacuum and the corresponding particle (more about this method can be found in [3, 10]) 

The LCQSR has wide range of applications to numerous problems of hadron physics. 
One of the promising ways for obtaining information about CKM matrix elements, as well 
as about wave functions, is studying the semileptonic decays. 

In this work we study — > (plu decay in LCQSR. This decay mode has been measured 
in experiments in [11]- [14]. Note that D — > <fi transition form factors are calculated in the 
framework of traditional 3-point QCD sum rules in [15, 16], but the results don't confirm 
each other. Therefore we decided to calculate D — > form factors using light cone sum 
rules as an alternative approach to the traditional sum rules. 

The paper is organized as follows. In section 2 we derive the sum rules for the transition 
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form factor. Section 3 is devoted to the numerical analysis and discussions and contains a 
summary of results and conclusions. 



2 Light cone sum rules for the D s —> <f> transition form 
factors 

We start by defining the form factors of D s — > <fi weak form factors in the following way 



((f)(P)\s^(l - 7 5 )c\D s (p Ds )) = -ie* fl {m Ds + m^)A 1 {q 2 ) 

+ -^®-wVf , (1) 

where q = po s — P is the momentum transfer, P and e are four momentum vector polar- 
ization of the vector meson, respectively, and po s is the four momentum of D s meson. 

In this section we derive sum rules for these form factors. In order to calculate the form 
factors of the semileptonic D s — > <p£u decay, we consider the following correlator function 

U lt (P,q)=iJ rf 4 xe^(0(P)T[s(x)7 M (l-75)c(x)c(O)(l-75)s(O)]|O) 

= r% - r + |^(2P + q), - r-^ + ir y e ^F . (2) 

The Lorentz invariant functions T°' ±,v can be calculated in QCD for large Euclidean p 2 D , 
to put it more correctly, when m 2 c — p 2 Ds 0, the correlation function (1) is dominated by 
the region of small x 2 and can be systematically expanded in powers of deviation from the 
light cone x 2 = 0. 

The main reason for choosing the chiral current c(l — 75)5 is that, in this case many 
of the twist-3 wave functions which are poorly known and cause the main uncertainties 
to the sum rules, can effectively be eliminated and provide results with less uncertainties. 
The chiral current approach has been applied to studying B — > n [17, 18], B — > rj [19] weak 
form factors. 

Let us discuss firstly the hadronic representation of the correlator. This can be done by 
inserting the complete set of intermediate states with the same quantum numbers of the 
current operator c(l — 75)5 in the correlation function. By isolating the pole term of the 
lowest pseudoscalar D s meson, we get the following representation of the correlator function 
from hadron side 

n (P „x _ (01^(1 -75)c|^)(^|c(l-75>|0) 
ll ^ q) ~ m* D .-(P + q y 

v (^(l- 75 ) C |fe)(fe|c(l-7 5 )a|0> 

Y rn 2 -(P + qy ■ [6) 
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For the invariant amplitudes T°' ±,v , one can write a general dispersion relation in the 
p 2 = (P + qf limit 



Y\q\ (P + q) 2 ) = J ds s _ P ^ S) +q)2 + subtr. , 



where the spectral densities corresponding Eq. (2) can easily be calculated. As an illustra- 
tion of this fact, we present the result for T 

p {0 \s) = ^^( mB + mv )A 1 (q 2 )5(s-m 2 Ds )+p^ h (s) . (4) 
m c + m s 

The first term in Eq. (4) represents the contribution of the ground state D s meson. In 
deriving Eq. (2), we have used 

(D s \c(l- l5 )s\V)= i - fl "' li 



m c + m s 



The second term in Eq. (4) corresponds to the spectral density of the higher resonances 
and continuum. The spectral density p^ h (s) can be approximated by invoking the quark 
hadron duality anzats 

p m {s) = p (o) Q cD {s _ So) 

So for the hadronic representation of the invariant amplitude we have 

r( o, = /».'< m B + m , ,~ p imco {s) suMr 

m c + m, m 2 n - (P + q) 2 J sn s - (P + q) 2 



Hadronic representations for other invariant amplitudes can be constructed in precisely the 
same manner. 

In order to obtain sum rules for the form factors A ± , A 2 , A and V, we must calculate 
the correlator from QCD side. This calculation can be performed by using the light cone 
OPE. The contributions to OPE can be obtained by contracting the quark fields to a full 
c-quark propagator, i.e., 

n M (P,g) = i I rf 4 xe^(0|s 7M (l -^)S c (x)(l - T5 )s(0)|0> 

= ij rf 4 xe^[Tr 7M (l - K)S c (x)(l - 75)^] (0|sFs|O) , (6) 

where T l is the full set of the Dirac matrices T l = (I, 75, 7a , 7q7 5, Pap), and 

1 }i + m c 



iS c (x) = iSi°\x)-ig s I -^e~ ikx f ' du 



(2k) 

H 9 l ux^G^iux)^ 

m 2 c — k 2 



2 (m 2 c - k 2 ) 



2 G^(ux)a^ 



(7) 
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Here, is the gluonic field strength, g s is the strong coupling constant and represents 
a free c-quark propagator 

From Eqs. (6)-(8) we see that, in order to calculate the theoretical part of the correlator, 
the matrix elements of the nonlocal operators between vector meson and vacuum states 
are needed. We see from Eq. (6) that, the contribution to the correlator comes only from 
the wave functions that contain odd-number 7 matrices. 

Up to twist-4, the meson wave functions containing odd-number of 7 matrices, and 
appearing in our calculations are: 



(<j>(P,X)\s(x)^s(0)\0) = 



P ^ Jo du ^ K " 2) + ^ 



l*^™ljy^C(u,^), (9) 



MP,A)|5(s)7 M 7fc*(0)|0) = U U - 2 J^) m ^elP a xp f due iuPx g { f> W) , 

4 y Tflfj) J JO 

(10) 

(<f)(P, X)\s(x)gG llu (ux)i-f a s(0)\0} = f<pm^p a (p u e x ± ^ - p^e x ±u )V(u,px) 

e x x 
px 

3 

(px)' 
e~x 



+ U m \—^^L-Vvgi^{u,px) 



+ U m $TZ^iVo l {ViiXv -PuX^{u,px) , (11) 



(<f)(P,X)\s(x)gG lxl/ (ux)i-f a -f 5 s(0)\0) = f^m^p^e^- p^e ±v )V(u,px) 



px 

In all expressions, we have used 

A _ e x x ( m\ 



+ f^rnl— Pa{PiJ,Xv - p u x^{u,px) . (12) 



px 2(px) X ^) +e± ^ 

gi» = 9l*> ~ —{Pv,Xv + PuX^) , (13) 
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where $y is the leading twist-2 wave function, while g±, g± , V are twist-3 and all the 
remaining ones are twist-4 wave functions. The notation used in Eqs. (11)— (14) is the 
following 

K(u, Px) = J Vae tPx(ai+ua3) K(a) , (14) 

where 



T>a = daida2das5(l — oc\ — ai — olj) . 



Inserting Eqs. (7) and (8) into Eq. (6) and using the definitions of <fi meson wave functions, 
the invariant structures r 0,± ' y take the following form 



/ 



2f (j) m (t> m c g ( f\u) 
Ax 

j f^m^m, 



I 



du 

du < 
Va 



(15) 



A? 



m 2 c miA(u) + 4(Pg)A 1 ${, i) (w) 



A 2 



A? 
A 2 



= J du ^-Um 4 ,m 



V V 

2 + 2, 



m 2 c m^A(u) + 4(Pq)A 1 $\ >(u 



f^m\m c uC{u) 
Af 



h m\m c {2 - u)C(u) 
Af 



(16) 



' 2$ - 2$ + ^ - ^ + - 

2 2 



where 



A? 



Ax = m 2 - (q + Puf , 

A 2 = m 2 c - [q + («i + wa 3 )P] 2 , 



(17) 
(18) 



Equating expressions of the invariant structures T 0,± ' v coming from QCD and phe- 
nomenological parts of the correlation function and making the Borel transformation with 
respect to (P + q) 2 in both parts, in order to suppress the contributions of higher states 
and continuum and also to eliminate the subtraction terms in the dispersion integral, we 
get the following sum rules for the D — > <p transition form factors: 



M<f) = 



i 



m c + m s 



fD 3 m 2 Ds m Ds + m^ 



,m 2 Ds /M 2 



2f (f) m 4> m c du 



9±{u) 



-s(u)/M 2 



U 



(19) 



x Y^m^m 
- 1 , 1 



1 1 



2 m l m c L -. A ( u ) 



-s(u) 



)/M 2 



uCW(u) (u)/M 2 

6 U 2 M 2 



rl 1 
brrArric I duVa8(so — s(k)) 



k 2 M 2 



2$(a) - 2$(a) + (a) - #(a) 



V V 

2 + 2 



-s{k)/M 2 



(20) 



The form factor A 3 (g 2 ) can be obtained from the exact result 



A 3 (q 2 ) = 



-Mi 2 ) 



2m^ ' 2m ( p 

and A (q 2 ) can be calculated from the following sum rule 



(21) 



A 3 (q 2 ) - A (q 2 ) 

X { W, 
,i i 



2 m c + m s q 2 1 * jM i 



f D ,m 2 Ds 2m - m 2 - g 2 
1 2 2 Z" 1 !., v 1 



3 -s(«) 



J<5 W 2 11 



f i / 2 2 2 2 

/ du —{m c — m^u —q 

J 8 U 



u - ■ u 2 M 2 



s(u)/M 2 _2 du 



(2-")g (il (") .-,(u)/^ 



1 



+ Um^rric duVa9(s - s ( k ))j^p 

s(k)/M 2 1 



W 2 M 2 

2$(a) - 2$(a) + - (a) 



V V 

-2 + 2 



(22) 



(m Ds + m<p) (m c + m s ) m i /M 2 
2m, f t 



x 1- 



where M 2 is the Borel parameter and 



Um, m J s 1 dugt\u)-^e^ M2 



(23) 



s(t) = 



m 2 c — qH + mitt 



t = < 



t 

u, or 



k — OL\ + M«3 , 

t = 1 — t , 
1 



2m 2 



ra«A + 9 - so) + V (so -ml- q 2 ) 2 - Ami (q 2 - m 
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3 Numerical analysis 



In this section we present our numerical calculation of the form factors A 1} A 2 , A and V. 
As can easily be seen from the expressions of these form factors, the main input parameters 
are the meson wave functions, whose explicit forms are given in [20] and we use them 
in our study. The values of the other input parameters appearing in sum rules for form 
factors are: m Ds = 1.9686 GeV, m s = 0.14 GeV, m c = 1.3 GeV, f Ds = (0.214±0.038) GeV 
[3], = 1.02 GeV. The leptonic decay constant of meson, which is f<f, = 0.234 GeV, 
is extracted from the experimental result of the — > decay [21]. The threshold 

s = (6.5 ±0.5) GeV 2 is determined from the analysis of the three point function sum rules 
for f Ds (see for example [3]). 

With the above-mentioned input parameters, we now proceed by carrying out our nu- 
merical analysis. The first step, according to sum rules philosophy, is to look a working 
region of the auxiliary Borel parameter M 2 , where numerical results should be stable for a 
given threshold sq. The lower limit of M 2 is determined by the condition that the terms 
M~ 2n (n > 1) remains subdominant. The upper bound of M 2 is determined by requiring 
that the continuum and higher state contributions constitute maximum 30% of the total 
result. Our numerical analysis shows that both requirements are satisfied in the region 
3 GeV 2 < M 2 < 4.5 GeV 2 . We should note that LCQSR for the form factors are reliable 
at the region q 2 < 0.4 GeV 2 . Moreover, we analyse the M 2 dependencies of the form factors 
Ai, A 2 , A and V at q 2 = GeV 2 and q 2 = 0.2 GeV 2 , for three different values of the 
continuum threshold, namely, Sq = 6.0, 6.5 and s — 7.0 GeV 2 . Our analysis shows that 
the form factors are practically independent of the Borel mass when M 2 varies between 
3 GeV 2 and 4 GeV 2 . Variation of the form factors in relation to the continuum threshold 
is also very weak. The results for all form factors change about 5% at q 2 = 0. Our final 
results for the form factors at q 2 = and s = 6.5 GeV 2 , are 

Ai(0) = 0.65 ±0.15 , 

A 2 (0) = 0.85 ±0.20 , 

A o (0) = A 3 (0) = 0.56 ±0.1 , 

V(0) = 0.90 ±0.20 . (24) 

It should be noted that in the region q 2 > 0.4 GeV 2 the applicability of the light 
cone QCD sum rule is questionable. In order to extend our results to the full physical 
region, we look for a parametrization of the form factors in such a way that in the region 
< q 2 < 0.4 GeV 2 , the above-mentioned parametrization coincides with the light cone 
QCD sum rules prediction. The most convenient parametrization of the q 2 dependence of 
the form factors is given in terms of three parameters in the following form 

= Z ( 21 2 F fl h ,2, 2 ,2 ■ (25) 

1 - a Fi {q /m Ds ) + b Fi {q /m D J 

The values of the parameters Fj(0), a Fi and b Fi are listed in Table (1). 

We proceed by discussing the uncertainties related to the input parameters and wave 
functions. We note first that the radiative corrections to the leading twist-2 function, which 
is calculated in [20], is about ~ 10%. As has already been noted, the results depend weakly 
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F(0) 


dp 


bp 




0.65 


1.36 


-0.31 


1 A Z 


0.85 


4.5 


5.55 


A, 


0.56 


0.13 


-0.46 


V 


0.9 


2.82 


1.51 



Table 1: Parameters of the form factors given in Eq. (25), for the D s decay in a three- 
parameter fit. We take the central values of the form factors for F(0). 



on the continuum threshold sq and Borel parameter M 2 , and the uncertainty due to these 
parameters is about 5%-7% in the working region of M 2 . Moreover, the results are also 
quite weakly dependent on the vector meson decay constant and fj, which results in 
an uncertainty about 5%. Additional uncertainty coming from the Gegenbauer moments 
are about ~ 10%. Summing up all these above-mentioned errors, the overall uncertainty 
in the values of the form factors is of the order of 17%. 
In the experiments, the ratios 



n = 



V(0) 



and 



Ai(0) ' T ' 2 A 1 (0) ' 

are measured. In the present work, within the framework of the light cone QCD sum rules, 
we get r\ = 1.57 ± 0.36 and r 2 = 1.30 ± 0.34. In Table (2), we present a comparison of our 
results with the existing experimental data and 3-point sum rule (3PSR). 

Using the parametrization of the D s — > <fi transition in terms of the form factors Ai, 
A 2 , V, A 3 — Aq, the differential decay width as a function of q 2 , in terms of the helicity 
amplitudes can be written as 



A 2 (0) 



dT 

dq 2 



G 2 F 



V r , 



m 



1927T 3 

dT L ^ (11 



H 2 + H 2 + + H 2 



+ 



dr. 



dq 2 dq 2 dq 



(26) 



where the indices in dTi/dq 2 and Hi denote the polarization of the <fi meson, X(m 2 D , mi, q 2 ) 



(m 2 D ^ +mi- q 2 ) 2 - Am 2 D mi, and 



H± = [m Ds + m^A^) =f —± V(q ) 



H 



2m (t ,y/q z 



[m 



D s 



m Da + 

m l ~ Q 2 )( m D s + m^A^q 2 ) 



\{m 2 Da ,ml,q 2 ) 
m Da + m<p 



(27) 



(2? 



The differential decay rate when the final state meson is transversally polarized is 
determined to be 



dT 



T 



dr. 



dr_ 



dq 2 dq 2 dq 2 



(29) 
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ri 


r-2 


E791 


2.27 ±0.35 ±0.22 


1.57 ±0.25 ±0.19 


CLEO 


0.9 ±0.6 ±0.3 


1.4 ±0.5 ±0.3 


E687 


1.8 ±0.9 ±0.2 


1.1 ±0.8 ±0.1 


E653 


2.3t^ ± 0.4 


2.it{|:5 ±0.2 


Average 


1.92 ±0.32 


1.60 ±0.24 


3PSR [15] 


2.20 ±0.85 


1.16 ±0.46 


3PSR [16] 


2.16 ±0.38 


-1.08 ±0.17 


Our Results 


1.38 ±0.44 


1.31 ±0.43 



Table 2: Comparison of our results for T\ and r 2 with the experimental results and 3-point 
sum rule. 

Integrating the differential decay widths over q 2 in the region from q 2 = to (m Ds — m^) 2 , 
we obtain 

T L = (1.52;°$) x 10- 14 GeV , 
T T = (2.21~°:?°) x ID" 14 GeV , 

and for their ratio, we get 

^ = (0.69 ± 0.44) , 
which is in good agreement with the existing experimental data 

| pi J = 0.72 ±0.16 , [21] 

\ T J eX p 

Using the value of the total decay width T Ds = 1.34 x 10~ 12 GeV [21] of the D s meson, 
we get the following result for the branching ratio of the D s — > <J)$i/ decay 

B(D S ^ fiv) = (2.78-};™)% , 

which is consistent with the experimental result 

B(D S -> <ftv) exv = (2.0 ± 0.5)% . 
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In conclusion, we calculate the form factors for the D s — > transition, in framework 
of the light cone QCD sum rules. We compare our results for the form factors with the 
existing calculations based on 3-point sum rules. Following this analysis, we then estimate 
the ratios of these form factors and compare them with the current experimental data, as 
well as with the existing theoretical calculations. Finally, we study the ratio T L /T T of the 
decay widths when <fi meson is longitudinally and transversally polarized, and the branching 
ratio. Our calculations on the above-mentioned quantities confirm that they are consistent 
with the existing experimental data. 
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